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ABSTRACT 

We use galaxy surface brightness as prior information to improve photometric redshift 
estimation. We apply our template-based photo-z method to imaging data from the 
ground-based VVDS survey and the space-based GOODS field from HST, and use 
spectroscopic redshifts to test our photo-z 's for different galaxy types and redshifts. 
We find that the surface brightness prior eliminates a large fraction of outliers by lifting 
the degeneracy between the Lyman and 4000 Angstrom breaks. Bias and scatter are 
improved by about a factor of two with the prior in each redshift bin in the range 
0.4 < z < 1.3, for both the ground and space data. Ongoing and planned surveys 
from the ground and space will benefit provided that care is taken in measurements 
of galaxy sizes and in the application of the prior. We discuss the image quality and 
signal-to-noise requirements that enable the surface brightness prior to be successfully 
applied. 

Key words: galaxies: distances and redshift - galaxies: photometry - methods: data 
analysis - surveys 
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1 INTRODUCTION 

Photometric redshifts ("photo-z's") are a tool for obtaining 
redshift and type information for galaxies for which broad- 
band photometric colors are available rather than spectro- 
scopic data. Generally, this technique is used to study pop- 
ulations of galaxies when the observational cost of obtaining 
spectra for these galaxies would be prohibitive. With the ad- 
vent of wide field imaging surveys, photo-z estimation has 
become an indispensable part of cosmological surveys. With 
imaging in three to five optical bands, ongoing and planned 
surveys aim to get photo-z's over hundreds or thousands of 
square degrees of sky. With the redshift and sky positions, 
surveys such as DE^|, KiDlfl LSST0, and Pan-STARRS|3 
plan to use galaxy clusters, galaxy clustering and gravita- 
tional lensing as probes of dark energy and other cosmolog- 
ical issues. The techniques and accuracy of photo-z's have 
become an area of active study since they play a critical role 
in deriving cosmological information from imaging surveys. 

Recently, two main approaches have been used to 
obtain photo-z's: empirical fitting and template fitting. In 
the former, a neural network, polynomial function, or other 
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empirical relation is trained using a subsample of galaxies 
with known (spectro scopic) redshifts, and then applied 



to the larger sample dConnollv et af] Il995l; [Vanzella et al 



2004; Firt h et all 120 03: Brodwin et al. 2006; Collistc r et al 
boon lOvaizu et all 120081: [D'Abrusco et all 120071 ; iLi et al 



2006; Abdalla et al.l 



2007 



iBanerii et al.l 120081 ). In the 

first a library of theoretical 



template-fitting method 
or empirical spectral energy distributions ("SEDs") are 
generated, and then they are fit to the observed colors 
of galaxies, where t h e redshift is a parameter that is fit 
llAmouts et al.l | l999l: [Benite zl lioool: iBolzonella et al.ll2000l: 



Babbcdgc et al 



| 2006l: iFeldmann et al.l 120061; 



2004: lllbert et al.l l200rT iBrodwin et all 
Mobasher et al.l 12007 



Brodwin et al. 



2006 

Wittman et all 120071 ). 



| 2007| : 

boon 



a 

be 



Margoniner fe Wittmanl 
As in the training set approach, 
subsample of galaxies with spectroscopic redshifts can 
used to help calibrate this procedure. In this study 
we use a template-fitting me thod, which makes use of a 
library of galaxy templ ates dBruzual A. fe Charlotl Il993l : 
iBruzual fe Charlotl 12003). to predict the colors of galaxies 
through a series of optical and infrared passbands. In 
principle the idea is simple: given an observed galaxy, 
compute the redshift z that makes each template T match 
the observed colors most closely, and then choose the best 
fit (T, z) pair. 

The reality is more difficult. With only a limited num- 
ber of colors, it's often impossible to tell whether a set of 
observed colors better matches one template which is at high 
z, or another at low z. This is known as the "color-redshift 
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degeneracy," and occurs because a red galaxy template at 
low redshift and a blue galaxy template at high redshift can 
look the same in a given set of filters. For the same reason, 
accurately determining the rest-frame color of a galaxy (the 
"fc-correction" ) is difficult: intrinsic astrophysical variations 
and evolution broaden scatter in the color-z relation, and 
the effect increases at higher redshift. With this paper we 
address the question of how to break this degeneracy. The 
most widely used method uses empirically m easured appar - 
ent magnitude-redshift distributions (see e.g. Benitcz 2000); 
in this paper, we introduce a new and complementary ap- 
proach to breaking the color-redshift degeneracy using the 
surface brightness (SB), which is the luminosity of an object 
per unit surface area. We will show that the SB is able to 
provide a strong constraint on galaxy redshifts, which should 
remain intact even for faint samples, and that this will not 
be the case for the magnitude prior. In order to incorporate 
prior statistical information about the galaxies under study 
into the redshift estimate in a consis tent way , we u se the 
Bayesian approach first introduced bv lBenitea (|2000T ) . 

Throughout this paper we use the AB magnitude sys- 
tem. 



2 BAYESIAN PHOTOMETRIC REDSHIFTS 

Our goal is to find, for each galaxy, the 2-D posterior distri- 
bution P(z, T\C, O), where z is the redshift of the galaxy, 
T is the "template parameter," which is a discrete variable 
corresponding to the galaxy type in our template library, C 
is the vector of fluxes from the data, and O is a vector of 
observables independent of the fluxes C for each galaxy. O 
could be any set of observables, e.g. size, brightness, mor- 
phology, environment. In the current paper we consider I- 
band apparent magnitude and apparent surface brightness. 
Using Bayes' theorem and the definition of conditional prob- 
ability, we can write 



<*) = 



P(z,T\C,0) 



P{C,Q,z,T) _ P(C\Q,z,T)P(Q,z,T) 



P(C,0) 



P(0)P(C) 



,(1) 



because we make the approximation that O and C are inde- 
pendent. By the same token, we can write P(C|0,z,T) = 
P(C\z,T), and then since P(z,T, 0)/P(0) = P(z,T\Q), 



P(z,T|C,0) = 



P{C\z,T)P{z,T\Q) 
P(C) 



(2) 



This is the equation on which our work is based. The pos- 
terior distribution P(z,T|C,0) is given in terms of the 
likelihood function P(C\z,T) and the prior distribution 
P(z,T\0); the prior encompasses all the prior knowledge 
we wish to use about galaxy morphology, evolution, envi- 
ronment, brightness, etc. Once the posterior distribution is 
determined, all the quantities of interest can be calculated; 
ideally, any study that wishes to make use of photo-z infor- 
mation would use the full P(z, T\C, O) directly. For the sake 
of simplicity, we evaluate the performance of our estimator 
using the mode of P(z, TjC, O), which is the "best" redshift 



Z(P = Pr 



(3) 



We can also compute the marginalized redshift estimate 
given by 



/ dTdzzP(z,T\C,0) 
J dTdzP{z,T\C,0) ' 



(4) 



(5) 



and the variance 

, _ 2 _ fdTdzz*P(z,T\C,Q) 2 
{ ' { ' ~ JdTdzP(z,T\C,0) { ' ■ 

As we will show, we find that P(z,T|C) is typically multi- 
modal with one minimum in \ close to the spectroscopic 
redshift and one or more "false" minima at other redshifts. 
Using z p from Eq. <(3j as the redshift estimator means that, 
without a prior, the algorithm can often generate a large 
error by choosing the estimate from one of the false minima. 
Including a prior deweights those minima and reduces the 
impact of the outliers (often referred to as "catastrophic" 
outliers), greatly improving the accuracy of the z v estimate. 
Ideally, however, any study using photo-z information should 
use the full posterior distribution P(z,T|C, O) in its analy- 
sis. 



3 TEMPLATE FITTING 

In order to calculate the likelihood function P(C|z,T), we 
initially model the colors of galaxies as a function of red- 
shift. We use the simple stellar population models from 
the GISSEL (Galaxy Isochrone Synth esis Spectral Evolution 
Library) spectral synthesis pac kage (|Bruzual A. fc Charlotl 
Il993l ; iBruzual fc Charlotl I2003T ) to derive a set of spectral 
energy distributions with ages from 50 Myr to 15 Gyrs. We 
parameterize the galaxy type as a one dimensional sequence, 
T, that expresses the age of the spectrum with T = 9 for 
galaxies with ages of 50 Myr through to T = for galaxies 
of age 15 Gyr. All SEDs are convolved with a set of filter re- 
spon se functions t aken from the GOOD S l|Giavalisco et al.1 
|2004) and VVDS l|Le Fevre et al.ll2004bl N ) surveys to model 
the color as a function of redshift over the redshift inter- 
val < z < 4 (see Figs. [T]and[2}. All input spectra are 
corrected for attenuation due to the intergalactic medium 
l|Madaulll995T ). We make no account for evolution in the 
SEDs as a function of redshift nor do we train the spectral 
models to correct for uncertainties within the photometric 
zero-points of the observed data; instead we choose a largely 
equivalent, and simpler, procedure of fitting a correction to 
the bias on the calibration sample for use on the (separate) 
test sample. 

Given a model of galaxy colors as a function of redshift 
and spectral type, Gi(z,T), whose values are the flux in 
the i th band for a given template T and redshift z we then 
calculate, for each galaxy, 



X 2 (z, T) = J2(Ci ~ aG t (z, T)fla\ 



(6) 



where P(C\z,T) = exp(-x 2 (z, T)/2), d are the observed 
fluxes, <Ji the photometric uncertainties and a is a scale 
factor that is a free parameter for each galaxy. The scale 
factor, a, means that the template fitting procedure by it- 
self does not take into account the overall brightness of a 
galaxy, but only its colors. Since we do not have an accu- 
rate way to estimate the true error in our template fitting 
procedure, once we have calculated the full 2-D likelihood 
function P(C\z,T), we scale the errors on the measured C 
so that Xmin(z,T) = 1. Once we have determined the prior 
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Figure 1. Flux density vs wavelength for the filters used in the GOODS and VVDS surveys. The ACS filters used are F435W (Johnson 
B), F606W (Broad V), F775W (SDSS i), and F850LP (SDSS z). The VVDS filters are Mould U, B, V, R, and I filters, as well as 
ESO/WFI J and K s filters. The detector response function for the ACS is indicated on the lower plot. 



distribution P(z,T\0), we merely have to multiply the two 
distributions to get the desired posterior distribution. We 
address how the priors are generated in detail in the next 
section. 

For the current work we utilize the galaxy catalogs from 
the GOODS and VVDS redshift surveys. For the GOODS 
data we utilize the B, V, i, z photometric passbands and for 
the VVDS the U, B, V, R, I, J, K s passbands (see Fig. H). 



4 PRIORS 

The choice of prior P(z,T\0) is key to a successful 
template-based photo-z study. The typical failure mode for 
a template-based photo-z code comes in the form of a mul- 
timodal probability distribution, where there is one peak at 
or near the correct redshift, but at least one other false peak 
in (z, T)-space (see Figs. [3] and [4}. The astrophysical origin 
of these multiple peaks in (z, T)-space (and in fact the abil- 
ity for photometric redshifts to perform as well as they do) 
comes from the Lyman (1216 A) and Balmer (4000 A) breaks 
in galaxy spectra. The transition of these breaks through a 
set of filters as a function of redshift provides the change 
in galaxy color that all photometric redshift techniques key 
into. If our photometric observations are sufficient to mea- 
sure the galaxy flux for both breaks then the templates usu- 
ally have the power to determine the galaxy redshift unam- 
biguously, and there is no degeneracy. This is often not the 
case for optical surveys. For example, the GOODS survey 
niters cover the range of approx. 3700-10000 A, so the Ly- 
man break is only observable for galaxies with 3 < z < 6, 
and the 4000 A break should be observable for galaxies with 
z < 1.5. With only sufficient spectral range to identify a 
break in the spectrum, but not to classify which break it 
corresponds to, the template fitting procedure will produce 
multiple peaks in the P(z,T) distribution. The VVDS adds 



the infrared filters (J and K s bands) which should help iden- 
tify the 4000 A break out to higher redshifts and, therefore, 
disambiguate the spectral breaks; however, at the depth of 
the /-band selected and / < 24 magnitude limited sample 
the infrared observations do not have sufficient signal-to- 
noise to constrain the infrared colors. In the surveys under 
study, there are almost no non-QSO/AGN galaxies above 
z = 2 that have spectra which meet our quality criterion, so 
most of the information available to our template fits comes 
from the 4000 A break. 

In the following subsections we examine how two differ- 
ent methods can break this degeneracy: first, our approach 
using an empirically-measured SB(z) relation, and then a 
widely-used method using the apparent magnitude. 

4.1 Surface Brightness Prior 

4.1.1 Tolman Test 

In this section, we consider the performance of using SB to 
break the color-redshift degeneracy. The use of SB in con- 
straining photo-z's has been proposed a number of times , 
for example, in the "/^-PhotoZ" method (|Kurtz et alj|2007h . 
which is calibrated with the SDSS. The authors use one- 
color photometry and surface brightness to get photometric 
redshifts for very red galaxies with z < 0.8. While their 
method works well the reddest 10-20% of galaxies, the con- 
straints worsen for bluer types, and thus it is probably not 
su itable for a general p hoto-z survey. The SB was also used 
in lWrav fc GunrJ (|200&t l as part of a larger framework to im- 
prove photo-z's for a low redshift sample of SDSS galaxies. 

The surface brightness (SB) of a galaxy is defined (in 
units of mag/arcsec 2 ) as: 

SB = m + 2.5 log A, (7) 
where m is the /-band magnitude and A is the angular 
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Figure 2. Galaxy templates fGISSEL IBruzual A. fc Charlodll993l: IBruzual fc Charloj|2003h transmission vs A, overplotted and normal- 
ized to the same value at 10000 A. The templates on the bottom, shown in red, are the reddest and oldest galaxy populations, while 
those toward the top are bluer and younger; the ages range from 14 Gyr (T = 0) to 50 Myr (T = 9). Convolution with the survey filters 
flattens out the emission lines, while the 4000 A and Lyman breaks are still pronounced, which are the primary features of the galaxy 
spectra that constrain photo-z's. 



area. The SB increases as the galaxy gets brighter and 
smaller. If we let x( z ) be the comoving distance to red- 
shift z, and recall that m ~ 2.5 log [x 2 (z)(l + z) 2 ] and 
log A ~ log {x~ 2 (z)(l + z) 2 ], then neglecting fc-correction, 
one can plug in to Eq. and derive the evolution of the 
SB with redshift: 



SB(z) = 10 log (1 + z) + const. 



(8) 



Equivalently, in flux units we get SB(z) oc (1 + z)~ 4 . This 
is how we would expect the SB to behave in an expanding 
universe for passively evolving galaxies for which there is no 
magnitude or size evolution with z. An interesting point is 
that SB(z) has no explicit dependence on the cosmology; in 
particular it is independent of the luminosity or angular di- 
ameter distances, unlike the magnitude-redshift relation. In 
fact, the SB(a) relation has been used to test the hypothesis 
that the u niverse is expanding , which is known as the "Tol- 
man test" (|Tolman| [l930, 193-i lLubin fc Sandagdl200ll ). In- 
stead of using the expected surface brightness-redshift rela- 



tion to test whether the universe is expanding, we use SB(z) 
to constrain galaxy redshifts. 



Just as in iLubin fc Sandagel (|200ll ). our galaxy sam- 
ple does not universally follow the simple law given 
in Eq. ((8]), because populations of galaxies evolve dif- 
ferently with redshift . The galaxy template library (see 
iBruzual A. fc Chariot] 1 19931 : IBruzual fc Chariot! |2003l , and 
Fig. [2} enables us to assign spectral types to each galaxy. We 
can, therefore, attempt to derive the magnitude fc-correction 
for each galaxy using its best fit spectral template and es- 
timate of the redshift. Fig. [5] shows the SB(z) relation for 
each galaxy type in our datasets, using the spectroscopic 
redshift information for the fc-correction. We can see that, 
in the panels with fc-correction, the redder galaxies are evolv- 
ing more passively and hence closer to the (1 + z)~ 4 slope, 
while the bluer galaxies are actively evolving in such a way 
that their SB doesn't change substantially with redshift; 
from this we infer that there must be significant evolu- 
tion effects that cause intermediate and blue galaxies at 
higher redshift to be brighter and/or smaller than passive 
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Figure 3. The two panels contain 1— (cyan), 2— (blue), and 3 — a (green) probability contours for the two galaxies in Fig. [4] Each 
galaxy has a multiply-peaked P(z,T) distribution where the template fitting approach fails to select the correct peak without a prior; 
the red triangle and magenta asterisk marks denote the values of (z) and z(x = Xmin)i respectively. For each galaxy, the upper panel 
shows how the SB prior breaks the degeneracy by eliminating the inaccurate high-redshift maximum, moving both {2) and z m - ln closer 
to the measured spectroscopic redshift, denoted by the dashed line. 
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Figure 4. Flux vs wavelength for the two degenerate galaxies shown in Fig. [3] The black curve is the high-resolution spectrum of the 
galaxy, the dashed blue and solid red curves are the template fits at high and low redshift, respectively, and the green bars are the 
observations in the broad band filters. The Lyman and 4000 A breaks are shown by dashed lines. For each galaxy, a P(z,T) distribution 
is computed by fitting the templates to the observed galaxy colors. However, without prior information, it is nearly impossible to tell 
which break is present in the spectrum, and therefore, whether the galaxy is at high or low redshift. 



evolution would indicate. Although they are constrained to 
lower red shifts (z < 0. 8), w e can see similar behavior in 
Fig. 3 of iKurtz et all {2007), where the lower left panel, 
containing the 10% reddest galaxies, follows the expected 
SB oc (1 + z)~ 4 relation most closely, with the strength of 



the SB(z) correlation decreasing as bluer deciles are consid- 
ered. 

In the plots with fc-corrections, the measured slopes are 
steeper for red galaxies; however, we use the non-fe-corrected 
SB to determine P(z,T|SB). The reason for this is that the 
effect of fc-correction is to remove the magnitude evolution 
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Figure 5. Surface brightness (SB) as a function of spectroscopic redshift, from space and from the ground, with and without magnitude 
fc-correction derived from the spectroscopic redshift z s and best-fit template. The panels are labeled by galaxy type (see Fig. [2j and 
best fit slope of the SB-log(l + z) relation (see Eq. {0). When observed from space, the passively-evolving red galaxies are close to the 
(1 + z)~ 4 power law predicted by Eq. {SJ, shown by the cyan line in each panel. Because the blue galaxies evolve with z, they have a 
shallower slope; this effect becomes more dramatic when fc-corrections are considered. No attempt has been made to compensate for the 
larger seeing in VVDS. 



with z, but the templates already provide this information, 
so once we have a best-fit template we don't gain anything 
from explicitly fc-correctfng; in fact, it results in a larger scat- 
ter in SB(.z), because the uncertainties in the k-corrections 
derived from photometric redshifts increase the scatter in 
the SB-redshift relation (due to the "catastrophic errors"). 
The fc-corrected plots in Fig. [5] have a small scatter because 
they have used the spectroscopic redshift in order to deter- 
mine the fc-correction, i.e. they use extra information that 
the photo-z procedure doesn't have. 



4-1.2 Calibration 

We calibrate the SB prior via the following procedure: first, 
we bin our galaxies into types T; each galaxy in the calibra- 
tion sample has a spectroscopic redshift z 3 and a measured 
(non-fc-corrected) surface brightness SB bs (computed using 
Eq. (0 from the /-band magnitude and size measure). For 
each type, we then do a (2<r-clipped) linear fit 

y(z,T)=A(T)log(l + z) + B(T); (9) 



these fits are plotted in Fig. [5] Then, using the fit in Eq. ©, 
for every galaxy i we can now compute what its SB would be 



GOODS-S, N(SB,z = 0,T) 
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Figure 6. JV(SB, z = 0,T) with 2cr-clipped Gaussian fits for the GOODS-S and VVDS datasets. Using the fit from Eq.[9]we have shifted 
the galaxy population in each template bin to z = and plotted the resulting distribution. 
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Figure 7. Histograms of P(z,T) for GOODS-S and VVDS datasets. The solid curves are exponential fits: P(z,T) 
(z 2 /2q(T)) exp [— z/zq(T)}; the parameter zq(T) is fit separately for each type. 



if the galaxy were at some redshift z instead of the observed 
spectroscopic redshift z B : 

SBi(z, T) = SB Ii0bs - A(T) [log (1 + z a ,i) - log (1 + 2)] .(10) 

If we then use this equation to move all the galaxies to the 
same redshift z, we can make a histogram iV(SB|z,T) from 
the SBi; iV(SBjz = 0,T), along with Gaussian fits, is shown 
in Fig. [6] Once we have these fits, we can pick an arbitrary 
redshift zo, and using Eq. (jTTJ} , compute SBi(zo, T) for every 
galaxy. We can further see that 



<SB( 20 ,T)) = A(T)Iog(l + *o) 



(11) 



(SB(2o,T)> 



+(SB obs -A(T)log(l + z s )}, 
A(T)\og(l + z ) + B(T), (12) 



where we have used B(T), the maximum likelihood estima- 
tor, for (SB obs - A(T) log (1 + z s )). 

The Gaussian fits to the histograms in Fig. |S] give us 
the width of the distributions, cr(T), for each type. Now we 
can find P(SB|z, T), which we take to be 



Pt(SB\z,T) = 



exp 



(SB obs -(SB( 2 ,T)» 2 



2a 2 (T) 



(13) 
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T Isophotal cr(T) Ellipse a(T) 






0.37 


0.48 


1 


0.30 


0.45 


2 


0.47 


0.54 


3 


0.51 


0.65 
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0.67 


0.65 


5 


0.46 


0.81 
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0.33 


1.0 
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0.43 


0.55 


8 


0.43 


0.57 


9 


0.62 


0.96 



Table 1. Effect of size measure on SB measurement in the 
GOODS data. For each template type, the scatter (in arcseconds) 
of the SB(z) relations around the best fits are shown for different 
size measures, on the left using isophotal area, on the right using 
ellipse area. Using the ellipse instead of isophotal area increases 
the scatter, which is why the precision of the SB measurement 
suffers when using the ellipse area. 



substituting from Eq. (|lip . we get 
P(SB\z,T) = ' 



(SB obs - A(T) log (1 + z ) - B{T)f 



(14) 



exp 



2a 2 (T) 

Finally, we use Bayes' theorem to get the prior distribution 
P{z,T) 



P(z,T\SB) = P(SB\z,T)- 



(15) 



P(SB) ' 

where we have plotted P(z,T) in Fig. [7] and fit it with an 
exponential function 



P(z,T) = 



■ exp 



(16) 



zUT)^\ z (T), 

We find P(SB) for both the GOODS and VVDS samples is 
well-fit by a Gaussian distribution. 

4-1.3 Size Measures and Seeing 

For different public datasets, different size measures are 
available. In the GOODS data, the ellipse axes, half-light ra- 
dius, (filtered) isophotal area, Gaussian FWHM, and Kron 
aperture were all cataloged. For the VVDS data on the other 
hand, only the ellipse axes were available. For the HST data 
we found that the filtered isophotal area above the analysis 
threshold (SExtractor output parameter IS0AREAF_IMAGE) 
gave the best results for our purposes when measuring the 
angular area of the galaxies, giving the tightest SB(z) re- 
lation. In Table [1] we show that the scatter of SB(«) for 
the space-based data increases by a large amount for all 
galaxy types when we use the A_IMAGE and B_IMAGE ellipse 
axis parameters, instead of IS0AREAF_IMAGE, to measure the 
angular area. 

We quantify the impact of the algorithm used to mea- 
sure area in terms of the impact of a degradation in the 
image quality due to seeing. In Fig. [8] we apply a transfor- 
mation to the galaxy sizes, 



\/ r a + r sccing' 



(17) 



for the red galaxies (T — [0, 1,2]) and then recompute the 
RMS scatter in the SB(z) relation for those galaxies. We 



can see that in order to make optimal use of SB information 
for photometric redshifts, one of the primary issues is that 
one needs a size measure that is precise even in the presence 
of noise and seeing: Fig. [8] shows that adding 1" of seeing 
to the isophotal area measurement increases the scatter in 
SB(z) by about 10%, but if the ellipse area is used for the 
measurement, there is about a 80% increase in scatter for 
the same amount of seeing. 

We also see in Fig. and Table Q] that the best-fit slope 
of the red galaxies when measured from the ground is much 
shallower than the best-fit slope when measured from space. 
Recall that we expect the blue and intermediate galaxies 
to have a shallow slope because they are actively evolving 
with z, but the red galaxies should have a steep slope even 
in the presence of seeing. However, the ground-based ob- 
servations in Fig. [5] clearly have a minimum SB cutoff at 
SB ~ 26.5 mag/arcsec 2 , below which no sources are de- 
tected. This means that the increased scatter due to the 
seeing and bad size measure naturally leads to a shallower 
slope, because larger scatter means that the points in the 
SB population which would be measured below the cutoff 
are instead not observed, while the points that scatter to 
larger SB are unaffected. The net result is similar to that 
which would be produced by a Malmquist bias, which will 
also affect any flux- limited SB(jz) measurement (our sample 
should be largely free from the Malmquist magnitude bias 
itself because, as illustrated in the right panel of Fig. [5] the 
spectroscopic sample of galaxies which we are using is well 
above the flux limit of the photometric data). 

Fig. O Table [T] and Fig. [S] together show the perfor- 
mance that we would need from a ground-based size mea- 
surement if we intend it to be useful for photo-z's. For the 
space-based data with the isophotal size measure, the see- 
ing starts to degrade the scatter around r see ing ~ 1", so 
ideally the seeing from the ground should be less than this 
in order to obtain the most constraining power. Regarding 
the sensitivity, we note that purely for the purposes of con- 
straining photometric redshifts, as long as the calibration 
sample is representative of the total sample, the photo-z 
measurement should be unbiased, regardless of what the SB 
detection threshold is. If for other purposes we want to mea- 
sure the slope of SB(.?) in an unbiased way, i.e. with a slope 
that is not made shallower by the combination of seeing and 
noise, a planned survey with a similar depth (z mc d ~ 0.7) 
from the ground should be sensitive enough so that the de- 
tection threshold for galaxies is below 26.5 mag/arcsec 2 ; this 
requirement will change for a deeper survey, but at the same 
time the SB detection threshold naturally becomes more 
sensitive as signal-to-noise increases. 

4.2 Magnitude Prior 

Another possible approach to solving the problem illustrated 
in Fig. |4] is to use the overall brightness of a galaxy as 
measured in one of the photometric bands; recall that in 
the template fitting procedure, in Eq. (|6]), only the col- 
ors and not the brightness of the galaxy are used. If we 
take O = I (the J-band magnitude) in Eq. p). then we 
can use this information to try to break the degeneracies. 
Sim ilar techniques have been used in previous studies (see 
e.g. iBenitezI l2000l: lllbert et all 120061: jBrodwin et all 120061 : 
iFeldmann eta! . 2006; Mobasher et allhoOTl among others). 
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In iBenitej <|2000l ) the author demonstrates how to self- 
consistently compute a magnitude prior of this type from 
the galaxy data under study. We do not take that ap- 
proach, and further make the simplifying assumption that 
P(z,T\0) — P(z\m), i.e. galaxy type is statistically inde- 
pendent of magnitude. We take our prior information about 
the magnitude distribution of ga laxies with redshift from the 
DEEP 2 survey (|Coil et al.ll2004 ). In that survey they found 



the following probability density function (PDF) for their 
spectroscopic sample: 

P(z\I) oc dN/dz oc z 2 exp(-z/z ), (18) 

where zo is a function of the J-ban d magnitude. Th e param- 
eter zo is fitted for in Table 3 of (Coil et al. 20041) and we 
extrapolate linearly to deeper redshifts from their fit. From 
the shape of P(z\I) distribution plotted in Fig. 1101 we can 
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anticipate that the effect of applying a magnitude prior will 
be to always favor the lower redshift minimum in a P(z, T) 
probability distribution. Since most galaxies really are at low 
redshift, we will see that this approach works well; however, 
using a magnitude prior exclusively would not be expected 
to work well if one is interested primarily in faint galaxies, 
since a magnitude prior does not provide much information 
about faint galaxies, where the prior has a long tail out to 
high redshift. We can see that the SB prior, plotted in the 
right panel, does not suffer from this effect, providing a com- 
plementary constraint for faint galaxy samples. 



5 RESULTS 

In this section we will analyze the results of applying our 
photo-z algorithm to each of these datasets with SB priors. 
We will also examine the bias and scatter in our measure- 
ment errors. 

5.1 Datasets 

We use the Great Observatories Origins Deep Survey 
South (GOODS-S) and the VIM OS VLT Deep Survey 
(VVDS) as our calib r ation samples jGiavalisco et al.1 12004| ; 
IVanzella et ail 120051 ; iLe Fevre et all l2004bl lal. l2005h . The 
ESO CDF-S master spectroscopic catalog contains 1115 
galaxies with spectroscopic redshifts, of which we use 603 
that have spectroscopic redshift confidence ^ 95%, i.e. we 
select those that have quality factors greater than or equal 
to "B," "3," and "2.0" from the following sources: VLT / 
FORS2 spectroscopy Vers ion 1.0; VIMOS VLT Deep Survey 
(VVDS) Version 1.0; and lSzokolv et al.1 (|2004h . The VVDS 
spectroscopic catalog contains 8981 objects from the VVDS- 
F02 Deep field. Of these, 4180 objects meet our quality cri- 



terion 95% confidence in the spectroscopic redshift), cor- 
responding to quality flags 3, 4, 9, 13, 14, and 19. We use 
2090 spectra to calibrate our method on the ground based 
data. We then take the GOODS-N field, which contains 1814 
spectra, and the other half of the VVDS data, to test our 
photo-z method on data which are distinct from the calibra- 
tion sample. All results testing the performance of the SB 
and combination SB+mag priors are based on these inde- 
pendent calibration and test data sets. 

The magnitude and size distributions of the sample of 
galaxies under study are shown in Fig. [9] We can see that 
while the spectroscopic samples are very similar, the total 
size and magnitude distributions of the galaxies are substan- 
tially different for the ground and space data. This is because 
the HST data that we evaluate have a substantially higher 
signal-to- noise ratio for the sizes and magnitudes. This sig- 
nificantly affects the measurement of SB in each sample as 
well as the performance of the SB and magnitude priors. 

In order to evaluate the performance of the template- 
fitting photo-z algorithm, and any improvements that ap- 
plying priors produces, we define as a figure of merit 

A=^£l, (19) 

1 + Zs 

which is the fractional error in 1 + z v , where z v is the pho- 
tometric redshift estimate from Eq. ((3J . For a population of 
galaxies, (A) and a a measure the photo-z bias and the RMS 
size of the fractional error in 1 + z. 

Instead of applying a correction to the photometric 
zero-points and t raining the templa tes using the calibration 
sample, as in e.g. lllbert et al.1 (|2006T ). we adopt a simpler and 
largely equivalent procedure of fitting the bias (ignoring out- 
liers) on the calibration set. First we find a polynomial fit 
for the bias F(z s ) = (z p — z s ) for all galaxies with jA| < 0.5 
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in each redshift bin in the calibration set, and then for every 
redshift in the test sample we apply the correction 



(20) 



We note that a correction of this type is only expected 
to work well in regions adequately sampled by the spec- 
troscopic redshift sample, while a correction to the magni- 
tude zero-points may be effective over a larger range. The 
black points and histograms in Fig. Illl and Tab. show the 
characteristics of the unweighted photo-z template fitting 
algorithm on the test sample after the correction has been 
applied. 



5.2 Performance without Priors 

The performance of the template-based photometric red- 
shift estimator without priors on the GOODS and VVDS 
data sets is comparable with that foun d by other groups 
|Mobasher et alj|2007l : lllbert et alj|2006j ). The dispersion in 
the photometric redshift relation, defined at the 68 th per- 
centile of the distribution, has a scatter of aes 0.08 for 
redshifts z < 1.3, and ass =Ss 0.14 at redshifts z > 1.3. If 
all of the data, including outliers, is considered, the disper- 
sion of the photo-z relation is substantially larger (> 0.45), 
which demonstrates the non- Gaussian nature of the photo- 
metric redshift errors, which arises due to the presence of 
the outliers or degenerate points (see Fig. Illfl . As noted by 
iMandelbaum et alj l|2007l ). the small bias in the uncorrected 
(i.e. no correction of the form Eq. (|20p ) photo-z's without 
priors is due to a fortuitous cancellation: the outliers which 
scatter to high z cancel the bulk of the distribution which 
is biased low. 

As discussed in Sec. [4] the dominant failure mode for 
this kind of photo-z procedure is mistaking a 4000 A break 
for a Lyman break, which is visible on the scatter plots in 
Fig. [TT] as a region in which z s < 1 but z v > 2; this out- 
lier population is composed of galaxies whose probability 
distribution P(z,T) has multiple peaks (e.g. the galaxies in 
Fig. [3J) . Another type of failure can happen when points 
on both sides of the rest-frame 4000 A break are not ob- 
served in the survey filters. We can see this in the VVDS 
where the performance of the unweighted algorithm suffers 
at redshifts at z > 1.2 due to the lack of z-band data: at 
these redshifts, the VVDS outlier fraction and aes increase 
sharply as the 4000 A break redshifts out of the /-band; 
unfortunately, the J and K s bands are much shallower than 
their optical counterparts and so do not help constrain high- 
z galaxy photo-z's. A similar phenomenon is visible in both 
VVDS and GOODS for galaxies with z < 0.4, due to the 
lack of deep it-band imaging. 



5.3 SB Prior Performance 

Fig. [12] shows a breakdown of the impact of the SB prior 
on the photo-z error distribution in different redshift and 
type bins; we can see that the SB prior improves bias and 
scatter on the independent test sample in almost all redshift 
and template bins. We characterize the scatter in the pho- 
tometric redshifts as the dispersion given at the 95 th and 
68 th percentiles of the distribution, which are denoted by 
solid ((795/2) and dotted (aes) lines respectively; the two 



measures of scatter differ in how sensitive they are to out- 
liers. After the SB prior is applied, we obtain (795/2 < 5% 
for the GOODS-N and ct 95 /2 < 8% for the VVDS data in 
each bin in the redshift range 0.4 < z < 1.3. The largest im- 
provement in (795 is in the low redshift bins z < 0.6 for the 
GOODS and z < 1.3 for the VVDS data, because the SB 
prior eliminates the largest number of outliers in these bins; 
we can see that outliers increase the scatter at z < 0.6 by a 
factor of 2-3 for the GOODS data and by a factor of 1.5-2 
for redshifts z < 1.3 in the VVDS data. At higher redshifts 
the photo-z algorithm performs better on the space-based 
data from GOODS, due to the higher signal-to-noise, better 
size measurement, and z band data. For the templates, the 
largest improvements in accuracy and precision are gained 
for the reddest (types and 1) and bluest (types 8 and 9) 
galaxies. For the very old and red galaxy types, this is to be 
expected, because they evolve passively with redshift. As we 
see in Fig. [S] this means that galaxy types that have steep 
slopes close to 10 log (l + z) provide the strongest redshift 
constraints. For the very young and blue galaxies, we see 
from Fig. [U] that the absolute SB distribution is very tight 
for these types, so even though they have a shallow SB(z) 
slope the prior is able to provide a constraint on their red- 
shift. 

We can gain some intuition about the physical reasons 
for the difference between the SB(z) relation for the red and 
blue galaxy types by considering the surface brightness of 
a galaxy population with characteristic luminosity Lo and 
size Rq: 



SB 



Lo (l + z) 



Rl (1 + z)2p (l + z) 4 ' 



(21) 



where the exponents q and p describe the luminosity and size 
evolution of the sample, respec tively. Studies of starburst- 
ing galaxies dDahlen et al.ll2007l ) have indicat ed that the size 
evolut ion of these galaxies is fit by p « — 1. IWillmer et al.l 
(2006) find that M* B for blue galaxies in the DEEP2 sur- 
vey brightens by approximately 1.3 per unit redshift, i.e. 
q ~ 0.5. Therefore, ignoring selection effects, we see that for 
blue galaxies a combination of size and luminosity evolution 
should partially cancel the geometric factor of (l + z) -4 . We 
observe this in our data: the general trend in Fig. [5] is that 
for the bluer galaxies, evolution and selection together pro- 
duce an effect that cancels the (1 + z)~ 4 evolution, while for 
the red galaxies, the effect is much less pronounced. This is 
consistent with our knowledge that red galaxies evolve much 
less than blue galaxies do. 

We can also see how the bias (A) is improved in each 
redshift and template bin. By eliminating the high-z out- 
liers, we achieve bias of less than 5% in each redshift bin in 
the range 0.4 < z < 1.3 for the GOODS data, and less than 
4% for the VVDS in the same range. When the SB prior 
is applied, the photo-z estimator continues to perform well 
outside this redshift range on the GOODS sample, while for 
the VVDS, there is larger bias at very low and very high 
z. This is due primarily to the fundamental limitations of 
the VVDS survey bands (no deep infrared data) we have 
discussed previously, as well as the issues that can affect the 
SB prior in ground-based observations: seeing, sensitivity, 
and size measurement (see Sec. 14 . 1 . 3 p . 

Table [2] shows the global bias ((A)), scatter (cta, (795, 
(768), and fraction of outliers (N(\ A| > 0.2)) for the SB, mag- 
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nitude, and combination priors for our GOODS and VVDS 
test samples. Both the magnitude and SB priors have com- 
parable effects on the fraction of outliers in the photo-z re- 
lation, reducing it by a factor of 2 before any further cuts. 
If the goal of an analysis is as small a scatter as possible 
then applying both the magnitude and SB prior produces 
the best results for both GOODS and VVDS. 

We also show in Table [2] a simple way to reduce the 



outlier fraction in the results of our photo-z algorithm. First 
we sort the galaxies on the value of the posterior proba- 
bility P(z,T\C,0), i.e. the final \ 2 , at the best-fit (z,T). 
Then the subsequent rows after the first in each subsec- 
tion of Table [5] are produced by cutting a certain fraction 
of the worst-fitting (highest \ 2 ) galaxies from the sample 
and then evaluating the bias and scatter of the remainder. 
For the GOODS data, this produces an improvement in the 
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scatter and outlier fraction as well as the bias; however for 
the VVDS, making this cut entails a tradeoff between bias 
and scatter. In both cases we can achieve an outlier fraction 
of < 2% if we are willing to sacrifice 20% of the sample. 

5.4 Applications to weak lensing 

Lensing tomography refers to the use of depth information 
in the source galaxie s to get t hree-dimensional information 
on the lensing mass By binning source galaxies 

in photo-z bins, the evolution of the lensing power spectrum 
can be measured. This greatly improves the sensitivity of 
lensing to dark energy in cosmological applications. The rel- 
ative shift in the amplitudes of the lensing spectra is sen- 
sitive to the properties of dark energy. It depends on both 
distances and the growth of structure, thus enabling tests of 
dark energy or modified gravity explanations for the cosmic 
acceleration. Any errors in the bin redshifts propagate to 
the inferred distances and growth factors and thus degrade 
the ability to discriminate cosmological models. 

The capability of lensing surveys to meet their scien- 
tific goals will depend on our ability to characterize the 
photo-z scatter, bias and fraction of outliers. Photometric 
redshifts must be calib rated with an appropriate sample o f 
spectroscopic redshifts ijHuterer et al]|2006l ; lMa et alJ l2005). 
This may be done more cheaply by using auto- and cross- 
corre lations of photometric and spectroscopic redshifts sam- 
ples ijSchneider et al.l l2006h . which can also be used to es- 
timate the redshift distribution for a galax y sample where 
the calibration data is incomplete (see also [Schneide r et al.l 
l2006h . 

For lensing power spectrum measurements, broad bins 
in photo-z are sufficient - bin width Az ~ 0.2 or larger. 
This makes photo-z scatter less of an issue, since wide area 
surveys average over million(s) of galaxies per bin. However 
residual bias in the estimated mean redshift per bin leads 
to an error in cosmological parameters; this systematic er- 
ror can dominate the error budget if the photo-z's are no t 
well characterized and calibrated. In Huterer et al. (2006): 
iMa et al.l (|2005l ^l ; [Newman et all (120061 ) it is shown that the 
next generation ground based surveys that cover ~ 1000 
square degrees require residual bias levels below 0.01, and 
our method is close to this goal in the 0.4 < z < 1.3 range; 
however the more ambitious surveys planned by LSST and 
SNAP require levels below 0.002. Such an exquisite control 
over photo-z biases requires improvements in photo-z tech- 
niques and very demanding spectroscopic calibration data. 

Our results on the surface brightness prior show that it 
will be valuable in eliminating outliers that can bias photo- 
z's. In the literature outlier clipping is often performed be- 
fore photo-z's are evaluated, but there is not a well estab- 
lished basis for how the clipping can be done with real data. 
This is part of the reason for why photo-z's do not per- 
form as well on data as expected from tests on simulated 
galaxy colors. With more realistic testing, we expect that 
the surface brightness prior will emerge as an essential part 
of photo-z measurements for lensing and other applications. 
Also of value is understanding the relation to galaxy type 
and redshift: for lensing it is permissible to use s ub-samples 
that h ave well-behaved redshift distributions (see lJain et alj 
l|2007l ) for an application of this idea to "color tomogra- 
phy")- 



Lensing conserves surface brightness, but it can intro- 
duce magnification bias in the apparent magnitude. Hence 
using magnitude priors can cause subtle biases in lensing 
measurements; e.g., behind galaxy clusters magnitudes are 
brighter than in the field. Thus photo-z's estimated with a 
magnitude prior would place galaxies that are behind clus- 
ters at lower redshift than galaxies that are in the field, 
and the effect would be more pronounced for more massive 
clusters compared to smaller ones, thus biasing cosmologi- 
cal inferences from lensing measurements. While the bias is 
expected to be small, at the percent level, it does argue in 
favor of surface brightness priors for precision lensing mea- 
surements. 

While we have discussed weak lensing induced shear 
correlations, similar conclusions apply to other lens- 
ing applications such as galaxy-shear cross-correlations 
(known as galaxy-galaxy lensing and discussed recently by 
iMandelbaum et al.l (|2007l )V For cluster lensing, i.e. making 
maps of the projected mass distribution using weak and 
strong lensing, the scatter in the photo-z relation is more 
important since only the source galaxies in the cluster field 
provide useful lensing information. 



6 DISCUSSION 

We have shown that using the surface brightness-redshift 
relation, we can improve template-based photo-z results for 
both ground and space data. Our main result, presented in 
Sec. shows how the SB prior is able to help eliminate 
much of the color-redshift degeneracy that is present for 
many galaxies in the GOODS and VVDS samples. When 
compared to a standard unweighted redshift estimator, SB 
priors reduce the fraction of outliers in the photo-z relation 
from 12-16% to ~6%. This results in a decrease in the scat- 
ter in the photo-z relation (as defined by the 95 th percentile 
of the distribution) by a factor of two or more, while the bias 
improves from a value of A ~ 0.1 to A ~ 0.05; almost all 
redshift and type bins improve in the VVDS data, whereas 
the GOODS data improves likewise for all types but with 
an especially large effect at redshifts z < 0.6. 

At redshifts z > 1.3 the limitations of the passbands 
used in the VVDS give rise to a bias that becomes more 
negative with increasing redshift. The utility of the SB prior 
is dependent on seeing and the size measures provided; in 
the VVDS data, the sole public size measure (ellipse area 
measured by SExtractor) is not precise enough to optimally 
measure the surface brightness-redshift relation. We expect 
that more sophisticated techniques, such as fitting a light 
profile, would reduce the sensitivity of the SB prior to seeing; 
in particular, we require a size estimator that minimizes the 
scatter in the SB(z) relation. If a ground-based survey can 
precisely measure the angular area of galaxies and achieve a 
seeing of 1" or less, with a representative calibration sample 
that has the same sensitivity, then such a ground-based sur- 
vey should be able to do almost as well as one from space 
with respect to constraining photo-z's using the SB. 

Since our current application utilizes spectroscopic data 
to calibrate and test the SB prior, we are necessarily re- 
stricted to studying bright galaxies. Such a regime natu- 
rally favors a magnitude prior; as Fig. \W\ shows, the bright 
galaxy population is almost exclusively at low redshift. As 



14 Hans F. Stabenau, Andrew Connolly, Bhuvnesh Jain 



VVDS 



GOODS-N 



Prior 


cut 


< A >±V# 


o~A 


1X68 


<xg 5 /2 


outl. 


JNo prior 


1 no 

JLUU 


1 1 Q _J_ 1 1 


4o.y 


f, 9 


117 
ll.i 


1^7 
10. / 




yy 


11 Q _l_ 1 1 

11. o it 1.1 


A fi 7 


( : ; 9 
O.Z 


i i c 
11.0 


1^7 
10. i 




95 


12.0 ± 1.1 


48.8 


6.3 


11.8 


16.1 




80 


13.5 ± 1.2 


50.0 


6.6 


12.7 


17.8 


A /f o rr 


1 nn 


n a 4- n ^ 


zo.u 


O. 1 


O.o 


7 1 
( . 1 






1 yl 4- n A 


OO 1 
ZU. 1 


£ -1 


o.o 


R R 
O.O 




95 


-1.7 ± 0.4 


17.2 


5.1 


5.0 


5.7 




80 


-3.3 ± 0.2 


10.0 


4.4 


3.8 


1.5 


SB 


100 


-1.4 ± 0.5 


24.0 


5.1 


5.7 


6.2 
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80 
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-2.4 ± 0.4 
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5.4 
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5.4 


5.2 
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95 
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3.2 
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-0.2 
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23.1 
14.1 
9.1 
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3.9 
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4.3 
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3.5 
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Table 2. Summary of photo-z performance results for unweighted, magnitude, SB, and combination priors for VVDS (left) and GOODS- 
N (right). The subsequent rows for each prior beyond the first show how cutting a certain fraction of the galaxies based on the combination 
goodness-of-fit of the templates and prior (i.e. the final posterior \ 2 value) enables one to trade improved scatter for slightly worse bias, 
eliminating a larger fraction of the outliers. The columns of the table are the fraction of sample remaining, (A) and uncertainty, a & , 
0£± at the 68 th percentile, at the 95 th percentile, and fraction of outliers, all expressed as percentages; we follow the convention and 
define an outlier as having |A| > 0.2. 



the galaxies become fainter, the magnitude prior will have 
less and less predictive power; unfortunately, since most 
galaxies are faint, this is precisely the galaxy population we 
are most interested in obtaining photometric redshifts for. 
We do not expect the SB prior to suffer from this effect as 
much; if the sizes of galaxies are measured carefully so as 
to minimize the scatter in the SB(z) relation, then the SB 
prior should retain its statistical power for faint galaxies as 
well. 

Based on these results we expect, therefore, that the use 
of SB in constraining photometric redshifts (whether using 
a template based approach, neural networks, or other em- 
pirical relations) can substantially improve the robustness of 
the application of photo-z's in cosmology. The details of the 
application and its impact will be application dependent. In 
some applications, su ch as WL tomography l|Huterer et al.l 
l200rj ; lMa etld1l2005l ), it is important for the measured red- 
shifts to be as unbiased as possible, while others, such as 
Baryon Oscillations, are more sensitive to the scatter in the 
measurement. Subsequent papers will address the question 
of optimal measures for size from ground-based data and the 
resulting improvements in the SB prior. 
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